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Использование пчелиных алгоритмов 
для решения комбинаторных задач 


Работа посвящена решению задачи разбиения графа. Рассмотрены 5 алгоритмов: итерационный, 
эволюционный, генетический, муравьиный и пчелиный. Разработана программа на ЭВМ, реализующая 
описанные модели биоинспирированных алгоритмов. Экспериментальные исследования доказали 
эффективность пчелиного алгоритма по сравнению с другими алгоритмами. 


Введение 


В работе рассматривается общий подход к решению некоторых комбинаторных 
задач на основе пчелиных алгоритмов, относящийся к классу биоинспирированных [1]. 
Данный класс алгоритмов разрабатывался в рамках научного направления, которое 
можно назвать «природные вычисления». В основе идеи этих алгоритмов лежит моде- 
лирование биологических процессов, роевого интеллекта [2]. 

В ходе проделанной работы были реализованы и исследованы следующие 
алгоритмы: 

Микроэволюционный [3]. 

Эволюционный [3]. 

Генетический [4]. 

Муравьиный [2]. 

Пчелиный [5]. 
работе проведены исследования и экспериментальные результаты сравнения 
работ предложенных алгоритмов на примере задачи разбиения графа. Поскольку в 
алгоритме не использовались специфические данные о задаче, то предложенный его 
можно применять и для других комбинаторных задач, в т.ч. и для задачи коммивояжера. 


м 


Постановка задачи разбиения графа 


Проблема разбиения графов возникает во многих научно-технических задачах, 
решение которых можно частично или полностью свести к решению на графовых 
структурах. К классу таких задач традиционно относят задачи математической физики 
на неструктурированных сетках, проблему балансирования загрузки распределенных 
вычислительных систем, проблему девиртуализации решения задачи на реальной парал- 
лельной ЭВМ и другие. Как правило, целью разбиения является получение незави- 
симых подмножеств вершин исходного графа с минимальным количеством связей 
между этими подмножествами. Формально это можно записать следующим образом: 

Мы имеем граф С = (Т.Е), где Г - это множество вершин графа, а Ё — 
множество его ребер. Количество вершин |И| равно и. Необходимо получить такое 
разбиение графа на р частей И,,Г.,...,Г», чтобы для всех пар 1, / = 1, 2,..., ри1*] 
выполнялись следующие условия: 

КОИ, =В 


ИСИ Ч...-Ы И, ЕГ. 
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Критерий — суммарное число связей между подграфами. Цель — минимизация 
критерия. Существуют различные приближенные методы решения этой задачи. Один 
из них — генетический алгоритм. 


Общие положения генетического алгоритма 


Генетические алгоритмы (ГА) — адаптивные методы поиска, которые в последнее 
время часто используются для решения задач функциональной оптимизации. Они ос- 
нованы на генетических процессах биологических организмов: биологические популя- 
ции развиваются в течение нескольких поколений, подчиняясь законам естественного 
отбора и по принципу «выживает наиболее приспособленный» (5агутуа| оЁ Фе ННез0), 
открытому Чарльзом Дарвином. Подражая этому процессу, генетические алгоритмы 
способны «развивать» решения реальных задач, если те соответствующим образом 
закодированы [6]. 

Основные принципы ГА были сформулированы Холландом (НоПапа, 1975), и 
хорошо описаны во многих работах [7]. В отличие от эволюции, происходящей в при- 
роде, ГА только моделируют те процессы в популяциях, которые являются существен- 
ными для развития. 

ГА используют прямую аналогию с механизмом, полностью имитирующим про- 
цессы, происходящие в живой природе. Они работают с совокупностью «особей» — 
популяцией, каждая из которых представляет возможное решение данной проблемы. 
Каждая особь оценивается мерой ее «приспособленности» согласно тому, насколько 
«хорошо» соответствующее ей решение задачи. Например, мерой приспособлен- 
ности могло бы быть отношение силы/веса для данного проекта моста. (В природе 
это эквивалентно оценке того, насколько эффективен организм при конкуренции за 
ресурсы). Наиболее приспособленные особи получают возможность «воспроизводить» 
потомство с помощью «перекрестного скрещивания» с другими особями популяции. 
Это приводит к появлению новых особей, которые сочетают в себе некоторые харак- 
теристики, наследуемые ими от родителей. Наименее приспособленные особи с мень- 
шей вероятностью смогут воспроизвести потомков так, что те свойства, которыми они 
обладали, будут постепенно исчезать из популяции в процессе эволюции [5]. 


Генетический алгоритм для решения задачи 
разбиения графа 


Определим алгоритм кодировки и декодировки хромосомы. Кодировка решения 
происходит следующим образом: имеется вектор длиной Н, где |Н| - количество вершин 
графа, Н;= 1,2,...,п. Тогда первые |Гл| ген характеризуют содержание первого подграфа, 
а следующие |Т)] ген — второго подграфа и т.д., как показано на рис. 1. 
аллельй _|1121|314151|161718|9 
ген Н; ОИ 9 


подграфы |; Г, Г; 


Рисунок 1 — Кодировка хромосомы 


Начальная популяция решений создается на основе стратегии «дробовика», т.е. 
случайной генерацией. Оператор мутации работает на основе случайного обмена пози- 
циями двух генов. Такая операция не нарушает условий. Опишем пример действия 
простого обмена. Исходное решение: 


ген НМ; |11|13121718151416 
аллель! | 1121314151617 |8 


584 «Искусственный интеллект» 32010 


Использование пчелиных алгоритмов для решения комбинаторных задач 6 К 


Решение после перестановки (переставлены значения 2 и 5 позиции): 


ген Н; |118121713151416 
аллельй | 1121341516718 


Во избежание нереальных решений используется упорядочивающий кроссинговер. 
Оператор редукции основан на элитном отборе — в популяции решений остаются хро- 
мосомы с лучшими целевыми функциями (ЦФ). Структурная схема работы ГА отобра- 
жена на рис. 2. 


Инициализация начальной популяции 


условие останова? 
да 
Окончание 


Рисунок 2 — Структурная схема работы ГА 


Временная сложность генетического алгоритма для данной задачи представляет 
2 
следующую зависимость — О(и”), где и — число вершин. 


Общие положения муравьиного алгоритма 


Основной идеей алгоритма является моделирование поведения муравьев. Колония 
представляет собой систему с простыми правилами автономного поведения особей. 
Однако, несмотря на примитивность поведения каждого отдельного муравья, поведение 
всей колонии оказывается достаточно разумным [8]. Таким образом, основой поведения 
муравьиной колонии служит низкоуровневое взаимодействие, благодаря которому, в 
целом, колония представляет собой разумную многоагентную систему. Взаимодействие 
определяется через специальное химическое вещество — феромон, откладываемое му- 
равьями на пройденном пути. При выборе направления движения муравей исходит не 
только из желания пройти кратчайший путь, но и из опыта других муравьев, инфор- 
мацию о котором получает непосредственно через уровень концентрации феромонов на 
каждом пути. Итак, концентрация феромонов определяет желание особи выбрать тот 
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или иной путь. Однако при таком подходе неизбежно попадание в локальный оптимум. 
Эта проблема решается благодаря испарению феромонов, которое является отрицатель- 
ной обратной связью [9], [10]. 


Муравьиный алгоритм для решения задачи 
разбиения графа 


Для начала определим свойства муравья. 

1. Каждый муравей обладает собственной «памятью», в которой будет хранить- 
ся список вершин ./;„, которые необходимо посетить муравью А, который находится в 
вершине 1. 

2. Муравьи обладают «зрением», прямо пропорциональным количеству внутрен- 
них связей в подграфе их, где Г — номер вершины, в котором находится муравей; / — 
номер подграфа, которому принадлежит вершина /. 

3. Каждый муравей способен улавливать след феромона, который будет влиять 
на вероятность выбора следующей вершины. Уровень феромона в момент времени { 
на ребре О, будет соответствовать т,(И). 

4. Вероятность перехода муравья из вершины 71 в вершину / подграфа у будет 
определяться следующим соотношением: 


О.Р 


Рук (0) а [2 
УБР ел, (1) 


ТЕ, 
Р‚ь (1) = 0, / Е Л 
где а, В — параметры, задающие веса следа феромона, коэффициенты эвристики. 


Выражение (1) является математической моделью «колеса рулетки». На рис. 2 ото- 
бражен пример «колеса рулетки» при выборе ребра из вершины 1. 


Эф — 


м Р(1,2)=28% 
ы Р(1,3)=15% 
ы Р(1,4)=33% 
м Р(1,5)=24% 


Рисунок 3 — Выбор ребра 


Количество феромонов, откладываемое муравьем на ребре (1, /), задается в сле- 
дующем виде: 


пб: 
а р’ =Р®, 


0, (1,7) ® Т, (1) 


где О — параметр, имеющий значение порядка длины оптимального пути, [4() — 
длина маршрута Т»(1). 
Испарение феромона определяется следующим выражением [11], [12]: 


2+0 =а- рух @+У Аг, (3) 


где т — количество муравьев, р — коэффициент испарения. 


(2) 
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В рассматриваемом муравьином алгоритме все агенты изначально распола- 
ь 3 
гаются в случайно выбранных вершинах. Временная сложность алгоритма — О(и”), 
где 1 — количество итераций, т — количество муравьев, и — число вершин. 


Пчелиный алгоритм 


Данный алгоритм моделирует поведение пчел в естественной среде. Идея пче- 
линого алгоритма заключается в том, что все пчёлы на каждом шаге будут выбирать как 
элитные участки для исследования, так и участки в окрестности элитных, что позволит, 
во-первых, разнообразить популяцию решений на последующих итерациях, во-вторых, 
увеличить вероятность обнаружения близких к оптимальным решений [5]. Приведем 
основные понятия пчелиного алгоритма: 

1. Источник нектара (цветок, участок). 

2. Фуражиры (рабочие пчелы). 

3. Пчелы-разведчики. 

Источник нектара характеризуется значимостью, определяемой различными пара- 
метрами. Фуражиры закреплены за источниками нектара. Количество всех пчел в этих 
участках больше, чем на остальных. Среднее количество разведчиков в рое составляет 
5 — 10%. Вернувшись в улей, пчелы обмениваются информацией посредством танцев на 
так называемой закрытой площадке для танцев [2]. 

Пусть решение представляет собой, как и в ГА, вектор Н. Областью поиска 
нектара для пчел будет являться пространство поиска решений, размерностью п! (коли- 
чество всех возможных перестановок вектора Н). Расположение источника нектара 
характеризуется конкретной перестановкой Н, решением. Таким образом, координатами 
источника является решение Н. Количество нектара на источнике обратно пропор- 
ционально ЦФ. Участок имеет размеры, где размер — количество решений, «близких» к Н. 
Близость между векторами определяется значением расстояния Хемминга между ними. 
К примеру, решения {5,2,7,3,4,1,6} и {5,4,7,3,2,1,6} являются «близкими», т.е. в про- 
странстве поиска они располагаются рядом, находятся на одном «участке». 

Приведём словесное описание алгоритма пчёл. 

1. Генерация участков для поиска нектара. 

2. Оценка полезности участков. 

3. Выбор участков для поиска в их окрестности. 

4. Отправка фуражиров. 

5. Поиск в окрестностях источников нектара. 

6. Отправка пчёл-разведчиков. 

7. Случайный поиск. 

8. Оценка полезности новых участков. 

9. Если условие останова не выполняется, то п. 2. 

10. Конец работы алгоритма. 

Таким образом, ключевой операцией алгоритма пчёл является совместное иссле- 
дование перспективных областей и их окрестностей. В конце работы алгоритма попу- 
ляция решений будет состоять из двух частей: пчёлы с лучшими значениями ЦФ элитных 
участков, а также группы рабочих пчёл со случайными значениями ЦФ. Зависимость 
временной сложности нчелиного алгоритма от числа вершин — О(и”). 


Экспериментальные исследования 


Экспериментальные исследования проводились на различных графах, размер- 
ность которых варьировала от 10 до 1000. Все рассматриваемые графы являются пол- 
ными связными. Для проведения исследований была создана программа на ЭВМ. 
Основные характеристики ЭВМ, на которой производились эксперименты: 
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— процессор — АМО Типоп, 1,6 ГГц; 

— ОЗУ - 768 Мб. 

В главном окне программы отображаются: значения лучших ЦФ итерационного, 
эволюционного, генетического и пчелиного алгоритмов, графики зависимостей ЦФ от 
времени, рисунок разбиения и внешние связи. На рис. 4 приведены полученные реше- 
ния различными алгоритмами. Здесь показано разбиение графа из 100 вершин на 5 под- 
графов, каждый из которых мощностью 20 вершин. Также отображены внешние связи 
между подграфами. На графике видно, что пчелиный алгоритм показал наилучший 
результат. В порядке уменьшения качества решений имеем следующий список: 

1. Пчелиный. 

2. Муравьиный. 

3. Генетический. 

4. Эволюционный. 

5. Итерационный. 

Экспериментальные результаты работ алгоритмов приведены в табл. 1. 


Таблица 1 — Результаты работ алгоритмов 


Алгоритм Целевая функция 
Итерационный 822 
Эволюционный 804 
Генетический 765 
Муравьиный 694 
Пчелиный 690 


Для каждого алгоритма в данном эксперименте условием останова является 10 с. 


Выводы 


В ходе проделанной работы была создана программа на ЭВМ, реализующая опи- 
санные алгоритмы. Экспериментальные исследования показали эффективность примене- 
ния алгоритмов роевого интеллекта — муравьиный и пчелиный. Предложенная кодировка 
решения для пчелиного алгоритма оказалась эффективной. На рис. 4 изображен пример 
работы алгоритмов для графа размерностью в 100 вершин. 


© Разбиение графа 
Граф Свойства алгоритмов Алгоритм Запуск Справка 
Число вершин графа 
ИИ 
Значение лучшей ЦФ (ИА) 
622 о 
Значение лучшей ЦФ (ЭА) 
ООО 
Значение лучшей ЦФ (ГА) 
765 

Значение лучшей ЦФ (МА) 
69 
Значение лучшей ЦФ (ПА) 
вириндиЩенО 


Запуск 


эо®эоо 


500 1000 1500 2000 2500 3000 3500 4000 4500 5000 5500 6000 6500 7000 7500 8000 8500 9000 9500 10000 
Время, мс 


[Пчелиный алгоритм 


Рисунок 4 — Тест с 100 вершинами 


588 «Искусственный интеллект» 372010 


Использование пчелиных алгоритмов для решения комбинаторных задач 6 К 


В программе показывается результат разбиения в виде рисунка, где в блоках — 
подграфы, содержащие вершины, а линии связи — внешние связи. Здесь также отобра- 
жаются графики работ алгоритмов. 

Для графов с размерностью до 100 вершин «хорошее» решение находится менее 
чем за 10 с. Метод пчелиных колоний лучше остальных продемонстрировал способно- 
сть выходить из локальных оптимумов за счет использования случайного поиска. С увели- 
чением числа вершин преимущество пчелиного алгоритма перед вышеперечисленными 
в плане качества решения увеличивается. При малых размерностях задачи (до 80 вер- 
шин) муравьиный алгоритм показывает лучшие решения. Это связано с применением 
эвристики. Применение эвристики дает отрицательный эффект для больших графов, так 
как зависимость временной сложности муравьиного алгоритма от числа вершин выражает- 
ся кубической зависимостью — и”. При решении задачи коммивояжера лучше метода 
пчелиных колоний выглядели модифицированный муравьиный алгоритм [13], [14] и ге- 
нетический алгоритм с модифицированным оператором кроссинговера [15]. 

Перспективным улучшением данного алгоритма является: 

1.Применение эвристик. 

2.Направленное движение пчел в поиске лучших источников нектара. 

3.Реализация нечеткого логического контроллера для динамического управления 
параметрами пчелиного алгоритма. 

4. Гибридизация алгоритма. 
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